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I present an effective classical theory which describes the 
non-perturbative real-time dynamics in hot gauge theories 
and has the potential for numerical implementation. 



I. INTRODUCTION 

At high temperature, the naive perturbation theory is 
well known to break down at soft momenta because of 
large collective effects. During the last decades, there 
has been a considerable effort toward the understanding 
of such effects and their resummation into effective per- 
turbative expansions, with a better infrared behaviour. 
(See Refs. |l,|| for some recent reviews and more refer- 
ences.) In particular, it has been recognized, by Braaten 
and Pisarski that the gauge-invariant resummation 
of the collective phenomena at the soft scale gT involves 
the whole series of "hard thermal loops" (HTL's): these 
are one-loop corrections to the amplitudes of the soft 
(k gT <C T) fields due to their scattering off the hard 
[k ~ T) plasma particles. (See also Refs. J| - 0.) 

As well known, however, infrared (IR) problems persist 
even after resumming the HTL's: in general, they are as- 
sociated with (quasi)static magnetic gluons, or photons, 
which are not screened by the collective phenomena at 
the scale gT. 

In Abelian theories, we encounter only mass-shell di- 
vergences, as, for example, in the calculation of the 
fermion damping rate ^| : the divergences occur because 
a charged particle which is nearly on-shcll can emit or ab- 
sorb an arbitrary number of soft, (quasi)static, magnetic 
photons (the latter are virtual quanta, so they actually 
describe collisions in the plasma) . The kinematics is very 
similar to the old "IR catastrophe" of zero-temperature 
QED and, as a matter of facts, the divergences can 
be indeed removed jl^,[ll| by an all-orders resummation 
of the soft photons, a la Bloch-Nordsieck |J . 

In QCD, the IR problems are more complicated, be- 
cause of the mutual interactions of the magnetic gluons[]. 



*Based on a talk given at TFT98, Regensburg, Germany, 
August 10-14, 1998 

1 For definitness, I shall mostly use a QCD-inspired terminol- 
ogy. But the present considerations apply to any hot Yang- 
Mills theory, in particular, to the electroweak plasma at high 
temperature (T > T c ). 



Perturbation theory breaks down at the scale g 2 T be- 
cause of IR divergences. Specifically, the usual connec- 
tion between the powers of the coupling constant and 
the number of loops is lost: starting with some given 
order in g (e.g., 0(g 6 ) in the case of the pressure p2[), 
one finds diagrams with any number of loops which con- 
tribute at the same order. In some cases, perturbation 
theory breaks down already at leading order: this is what 
happens for genuinely non-perturbative phenomena like 
the static magnetic screening |L2| or the anomalous viola- 
tion of the baryon number in the hot electroweak theory 
HQ (see also below). 

For static quantities like the pressure, there exists a 
well defined procedure for going beyond the perturbation 
theory: this is lattice QCD at finite temperature [15|, 
possibly supplemented by dimensional reduction jl6], 1 7 1 . 
But the standard lattice simulations, as formulated in 
imaginary time, are not well suited for computing real- 
time correlation functions, like the baryon number viola- 
tion alluded to before. It is my purpose in this paper to 
describe a semi-classical approximation which correctly 
reproduces the non-perturbative real-time dynamics to 
lowest order in g and has the potential for numerical im- 
plementation (a concise presentation of this method has 
been previously given in Ref. [[l8|). 



II. THE CLASSICAL APPROXIMATION 

Let me first explain why one expects the classical ap- 
proximation to be appropriate for the problem at hand. 

The basic idea is not new (see e.g. [[l3|): the non- 
perturbative phenomena are associated with long wave- 
length (A <; l/g 2 T) magnetic fields, which, because of 
the Bose-Einstein enhancement, have large occupation 
numbers: 



N (k) 



1 



1 



T 

I >>1 



for fc«T, (2.1) 



and should therefore exhibit a classical behaviour. To see 
this in a simple way, let us put back Planck's constant h 
in Eq. (2.1) and compute the average energy per mode 
in thermal equilibrium: 



e(k) 



hk 



1 



T 



as hk < T. 



(2.2) 



As h — > 0, we recover the classical equipartition theorem, 
as expected. But Eq. (|]^) also shows that the relevant 
inequality is hk <C T, so that the classical limit (h — > 



1 



at fixed k and T) is actually equivalent to the soft 
momentum limit (k — > at fixed h and T). 

This observation is useful since we know how to per- 
form real-time lattice simulations for a classical thermal 
field theory: All we have to do is to solve the classical 
equations of motion for given initial conditions, and then 
average over the classical phase space with the Boltz- 
mann weight exp(— (3H). Since the initial conditions (say 
0(x) and 0(x) for a scalar theory) depend only on the 
spatial coordinate x, the phase space integration is ac- 
tually a three-dimensional functional integral, which can 
be implemented on a lattice in the standard way Jisl l- 

On the other hand, the classical approximation is cer- 
tainly not appropriate at high momenta, k <^ T, and this 
sets limits on its applicability (see below). 

To be more specific, let me consider classical Yang- 
Mills theory at finite temperature. In the temporal gauge 
A® a = 0, the independent degrees of freedom are the vec- 
tor potentials A z a (x) and the electric fields ( "the canoni- 
cally conjugate momenta") E l a (x), with the Hamiltonian: 



H- 



YAI 



= ~ / d 3 x {E a • E Q + B a • B a } . 



(2.3) 



In terms of these, the classical equations of motion read: 

d A° = -Et, d E° = ),!!,,)". (2.4) 

together with Gauss' law which in this gauge must be 
imposed as a constraint: 



D • E = 0. 



(2.5) 



The thermal pha se s pace is defined by the initial con- 
ditions for Eqs. (2.4), which I denote with calligraphic 
letters: £f (x) and Af (x). Then, the canonical partition 
function reads: 



J YM 



(2.6) 



whe re TI ym is expressed in terms of the initial fields, cf. 



Eq. (2J3). The classical thermal correlation functions can 
be obtained from the following generating functional: 



[Jf] = J m?VA18{V-£) 



x exp < —fiH-YM + / d*xJ? (x)A?(x) \ , (2.7) 



where A l a (x) is the solution to Eqs. (2.4) with the initial 
conditions {£f,.4f}. For instance, the free (g = 0) two- 
point correlation is obtained as: 

D° 3 cl (k) = (6ij - k t k 3 )Po(ko,k)N cl (k ), 

p (k , k) ee 2^e(fc ) 5{k 2 - k 2 ), N cl {k Q ) ee L , (2.8) 

KQ 



which is simply the soft (or classical) limit of the cor- 
responding Wightman functions^]: D^(k) ~ D^(k) ~ 
(T/k )p (k) EE D^ l {k). Other thermal expectation val- 
ues are defined similarly. For instance, the study of 
the baryon number violation at high temperature fTs) ] 
requires the quantity ((AB(t)) 2 ) , where 



AB(t) oc 



dx I d 3 xF^Fr, 



(2.9) 



(with Fg v ee {l/2)e^ pX F^ x ) is proportional to the 
change in the Chern-Simons number. And actually there 
have been attempts to compute {{AB{t)) 2 ) via lattice 
simulations of the equations above [|l9). 

Recall, however, that the classical d escri ption isjus- 
tified only at soft momenta (cf. Eqs. (2d) and fl2.2|) ), 
which grea tly restri cts its applicability. In general, equa- 
tions like (2J3) or ( |2.7| ) will be afflicted with ultraviolet 
(UV) divergences which occur because of the replacement 
of the correct, quantum distribution N (k) — l/(e l3k — 1) 
with the classical one N ct (k) = T/k (cf. Eq. Re- 
call, for instance, the famous "ultraviolet catastrophe" 
of Raylcigh and Jeans: The classical estimate for the en- 
ergy density of the black body radiation, namely (A is an 
ad-hoc UV cutoff): 



Ed/V 



d 3 k 
J2nf 



s c i(k) = T 



d 3 k 
(2^)3 



oc TA 3 , (2.10) 



is obviously wrong (since UV divergent), in contrast to 
the quantum result: 



E/V 



d 3 k 



(2tt) 3 e? k - 1 
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(2.11) 



which is finite because the large momenta k 3> T are 
exponentially suppressed by the Bose-Einstein thermal 
distribution. 

Similarly, the Debye screening mass is linearly UV di- 
vergent when computed in the classical theory. Recall, 
indeed, the correct, quantum, expression for the Debye 
mass to leading order in g : 

m% = -^rdkk^ = ^l, (2.12) 



JO 



dk 



where Ca = N for SU(A). In the classical theory, this is 
replaced by: 



._>,,, _ 9 2 Ca f A dkk 2^N cl _ g 2 C A TA 



/(A) = -■ 



dk 



(2.13) 



where the UV cutoff A has been introduced as an upper 
momentum cutoff. 



2 Recall the definition of the Wightman functions Q 
D>(k) = p(k)[l + N (k )] and D < (k) = p{k)N {k ). 
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Thus, the classical approximation, as formulated 
above, makes sense only for quantities which (at least, 
to some approximation) are not sensitive to the hard 
(k ~ T) plasma modes. Such quantities exist indeed: the 
magnetic screening mass is an example. More generally, 
all the static, non-perturbative, correlation functions of 
the magnetic fields are correctly reproduced by the clas- 
sical Yang-Mills theory, to lowest order in g. (In higher 
orders, the quantum corrections to the thermal distribu- 
tion start to play a role.) 

The easiest way to see this is to note that, in its mag- 
netic sector, the classical Yang-Mills theory is the same 
as three-dimensional QCD, a theory which is well known 
to describe the leadin g IR behaviour at high-T (see, e.g., 
HOI). Indeed, Eq. (§J) can be rewritten as: 



J YM 



= / VAt VAf 



cx P j-^y d 3 x 



(2.14) 



where the Aq components of the gauge fields have been 
reintroduced as Lagrange multipliers to enforce Gauss' 
law. As anticipated, the magnetic sector in Eq. (2.14) 
is equivalent to three-dimensional QCD with dimcnsion- 
ful coupling constant 173 = g 2 T. Incidentally, this shows 
that the "dimensional reduction" at finite temperature 
(cf. and Refs. therein) is a reflection of the classi- 

cal character of the leading IR dynamics. Thus, the vari- 
ous lattice investigations of QCD3 in the literature (e.g., 
the calculation of the magnetic mass pQ] , or of the non- 
perturbative correction, of 0(g 6 ), to the pressure |21|) 
can be also interpretated as classical calculations. 

On the other hand, we shall shortly see that the clas- 
sical Yang-Mills theory is not appropriate for the time- 
dependent non-perturbative correlations, which remain 
sensitive to the hard plasma modes. Accordingly, one ex- 
pects fl4|,^2| the aforementioned classical studies of the 
baryon number violation [JL9| not to be reliable. 



III. THE ROLE OF LANDAU DAMPING 

The reason for such a failure is that the time-dependent 
magnetic fields are sensitive to the hard thermal modes, 
via Landau damping Q. To see this, consider the propa- 
gator of a soft magnetic gluon in the HTL approxima- 
tion (that is, after resumming the effects of the hard 
modes to leading order in g). This reads D^(w,k) = 
(Sij - kikj)A T (uj,k), wherefl: 



A r (w,fc) 



j 2 - k 2 -U T (ui,k) ' 



(3.1) 



The retarded prescription is implicit here: lo = lo + i0~* 



and Ht(lo, k) is the corresponding HTL. 

In the static limit, Hr(0, k) — 0, and there is no sen- 
sitivity to the hard modes. More generally, all the HTL 
corrections vanish in the static limit, except for the term 
yielding Debye screening 0J^]. 

For non-zero lo, IIt(w, fc) is parametrically of the order 
m 2 D (u>/k), which is of 0{g 2 T 2 ) if lo ~ k: that is, for 
large enough frequencies, the time-dependent magnetic 
fields are screened as efficiently as the electric ones, and 
decouple from the non-perturbative physics. 

As lo — + 0, the polarization tensor is dominated by 
its imaginary part, which describes the absorbtion of the 
soft gluon by a hard thermal particle (Landau damping): 



n 1 u 

n T (aj <C fc) ~ — 1 — m n — 
4 fc 



(3.2) 



Accordingly, the magnetic propagator at lo <C k reads: 



A t (lo < fc) 



i (ttlo /4fc) m 2 D ' 



(3.3) 



which strongly suggests that non-perturbative phenom- 
ena can be associated only with frequencies which are low 



enough: 



< 



k 3 /r 



For fc ~ g 2 T, this implies lo ^ g 4 T. 



That is, the non-perturbative physics involves very soft 
and quasistatic magnetic fields, with spatial momenta of 
0(g 2 T) and even softer frequencies, of 0(g 4 T) @|. This 
conclusion is further supported by the power counting 
analysis in Ref. [fTif . 

Note that the inclusion of the hard modes effects (in 
the form of the HTL) was essential for the previous 
analysis. In particular, the typical time scale for non- 



perturbative phenomena r 



rir 



,/fc 3 - l/(g 4 T) is pro- 



portional to the Debye mass [Q, a quant ity w hich is 
dominated by the hard modes (cf. Eq. ( [2.12 )). We 
conclude that the dynamics of the time-dependent, non- 
perturbative, magnetic fields is strongly influenced by 
their scattering off the hard, quantum, modes, and can- 
not be described within the classical Yang-Mills theory. 

But this does not mean that the classical approxima- 
tion by itself is totally useless: we just have to be careful 
to apply it to the soft fields alone. That is, we have 
to treat hard and soft modes separately ]23] j: The hard 
(fc ~ T) modes are quantum but perturbative, so they 
can be explicitly integrated over in perturbation theory, 
to obtain an effective theory for the soft (fc ^ gT) modes. 
The resulting effective theory is non-perturbative, but for 
it the classical approximation is valid, so it can be studied 
numerically on a classical, three-dimensional, lattice. 

An essential step in this strategy is to clearly distin- 
guish between hard and soft degrees of freedom (to avoid 
overcounting, and to provide an ultraviolet cutoff to the 
effective theory for the soft modes). Loosely speaking, 
this requires an intermediate scale A, with gT « A « T, 
which should act as an infrared cutoff for the hard modes 
and as an ultraviolet cutoff for the soft ones, and which 
should cancel in the calculation of physical quantities. 

To lowest order in g, the effective theory for the soft 
modes is essentially known: this is the HTL-resummcd 
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theory by Braaten and Pisarski j|] , which includes all the 
one-loop amplitudes with soft external lines and hard 
loop momenta. And actually it has been already pro- 
posed by Bodeker, McLerran and Smilga (23] to use 
the HTL effective theory as a classical theory for non- 
perturbative calculations in real time. However, the prac- 
tical implementation of this suggestion met with several 
technical difficulties: (i) There is no intermediate scale 
A in the original formulation of the HTL theory: in- 
deed, this was not necessary for perturbative calcula- 
tions, where one can avoid overcounting by adding and 
subtracting the HTL Lagrangian to, and respectively 
from, the tree- level Lagrangian (ii) The HTL La- 
grangian is non-local 0,^], so it is a priori not obvious 
how to construct the associated classical thermal field 
theory (cf. the discussion in Sec. II). 

Actually, problem (ii) above can be overcome by using 
one of the local formulations of the HTL theory |||| . In 
what follows, I shall use the formulation in terms of ki- 
netic equations, as developed in Refs. || (see also 0). 
This has the advantage to clearly emphasize the physical 
content of the HTL's (see also Q]). Besides, we shall see 
later that problem (i) above — namely, the implementa- 
tion of the intermediate scale A — can also be solved in 
the framework of the kinetic theory Eq] . 



IV. KINETIC EQUATIONS FOR HTL'S 

I consider a purely Yang-Mills plasma in thermal equi- 
librium at a temperature T which is high enough for the 
coupling constant to be small: g <C 1. The typical exci- 
tations are "hard" gluons, with momenta k ~ T. Such 
gluons can develop a collective behaviour over a typi- 
cal space-time scale ~ 1/ffT, which is large as compared 
to the mean interparticle distance ~ 1/T. This results 
in long- wavelength (A ~ ^ 1/T) colour oscillations 

which are most economically described in terms of kinetic 
equations ^|). In these equations, the hard (k ~ T) 
gluons are represented by their average colour density 
5N a (k, x) to which couple the soft (i.e., long- wavelength) 
colour fields A£(x). The relevant equations read: 



(D v F">) a (x) = 2gC A 



d 3 k 
(2Vp 



v" SN a {k,x), 



(»-D I Uffl'(k )I )^jv.E a (i) 



dN 
dk 



(4.1) 



where = d^ L + igA£T a is the covariant derivative, 
E l a = F*° is the chromoelectric field, and v** = (1, v), 
with v = k/fc denoting the velocity of the hard particles 
(k = |k|, so that |v| = 1). 

The Yang-Mills equation above involves the colour cur- 
rent induced by the hard particles: 



f a (x) = 2gC A 



d 3 k 
(2tt) 3 



5N a (k,x). 



(4.2) 



By solving the second Eq. (4.1) (which may be seen as a 
non-Abelian generalization of the familiar Vlasov equa- 
tion |l],[7|,m ) ! we can express this current in terms of the 
gauge fields A%, and thus obtain an effective Yang-Mills 
equation which involves the soft fields alone: 



D V F^ =mj D — 



dn 



An v ■ D 



(4.3) 



Here, the angular integral J dfl runs over the orientations 
of v, a nd m j is the Debye mass, Eq. (2.12). 

Eq. (4.3) describes the propagation of soft colour fields 
in the high-T plasma. The hard particles are not ex- 
plicit anymore, since they have been integrated to yield 
the induced current in the r.h.s. By expanding this cur- 
rent in powers of the gauge fields one generates || all 
the HTL's of Braaten and Pisarski. However, because of 
the non-local structure of this current ( note the covari- 
ant derivative in the denominator), Eq. (4.3) is not very 
convenient for the construction of the thermal partition 
function [f23| . 

To this ai m, it is preferable to go one ste p ba ck and 
replace Eq. (4.3) with the coupled system (4.1) of lo- 
cal equations. There is a price to be payed for that: 
in addition to t he gauge fields A£(x), the local descrip- 
tion in Eq. (4.1) also involves the average colour density 
5 N a (k, x) , which can be seen as an "auxiliary field" . Still, 
when working with a local theory, we are in a better po- 
sition to look for a Hamiltonian formulation, as I discuss 



now. 



V. THE CLASSICAL EFFECTIVE THEORY 

The first step is to recognize, on the second Eq. (f4.1|), 
that the v and fc-dependence can be factorized in 
5 N a (k, x) by writing: 

SN a (k, a:) = -gW a (x, v) (dA /dfc). (5.1) 

The new functions W a (x,v) satisfy the equation: 

(v ■ D x ) ab W b (x,v) = vE (a:), (5.2) 

which is independent of k since the hard particles move 
at the speed of light: |v| = 1. By using Eq. (jO]), the 
induced current can be written as: 



.//» = '"/» / ^v»w a (x,v), 



(5.3) 



where the radial integration has bee n exp licitly worked 
out to yield the Debye mass (cf. Eq. Q2.12D ), 

The Hamiltonian formulation of the effective theory 
pl[ involves the auxiliary fields W a (x,v) together with 
the soft gauge fields A%(x). [Note that a different Hamil- 
tonian formulation of the HTL theory has been previ- 
ously given by Nair in terms of some other auxiliary 
fields.] In the temporal gauge Aq = 0, the independent 
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degrees of freedom are Ef, Af and W a , and the c orre - 
s pon ding e qua tions of motion can be read from Eqs. ( LI ) 
(|5.2D and (5.3) above. Specifically: 



/A O 
— Vi W a (x,v), 

(c\, + v • D) ab W4 = v-E Q , (5.4) 
together with the constraint expressing Gauss' law: 

G a (x) = (D • E) a +m 2 D I — W a (x, v) = 0. (5.5) 

J 47T 



Note that Eqs. (5.4) are not in canonical form: this is al- 
ready obvious from the fact that we have an odd number 
of equations. Still, it can be verified that these equa- 
tions are conservative; the corresponding, conserved en- 
ergy functional has the gauge-invariant expression Qj: 



H = - 
2 



d d xi E, 



+ m 2 D 



47T 



B„ ■ B„ 



W a {x,w)W a {x,w) 



(5.6) 



Moreover, in the gauge Aq = 0, the functional (5.6) acts 



also as a Hamiltonian, that is, as a generator of the time 
evolution. To see this, we need the following generalized 
Poisson brackets, as first introduced by Nair ||: 

{£?(x),4(y)} =-5 ab M (3) (x-y), 

{Ef (x), W b (y, v)} = Vi 6 ab 6^{x - y) , 
ml{W a (x,v)AV b ( y y)} = (gt bc W c + 

+ (vD x )° 6 ) ^(x-y)5(v )V '). (5.7) 

Here, <5(v, v') is the delta function on the unit sphere: 

dn 



-in 



*(v,v')/(v) = /(v'), 



(5.8) 



and all the other Poisson brackets are assumed to vanish. 
We also assume standard properties for such brackets, 
namely antisymmetry, bilinearity and Leibniz identity. It 
is then easy to verify that (a) the Poisson brackets (5.7) 
satisfy the Jacobi identity (as necessary for consistency) 
and (b) the equations of motion (|5.4j) f ollow as canon- 
ical equations for the Hamiltonian ( |5.6| ). For instance, 
d W a = {H, W a }, and similarly for Ef and A^ . 

Note also that the Hamiltonian in Eq. (^1]) is remark- 
ably simple: the piece involving the auxiliary fields Wq 
is simply quadratic ! In fact, all the dynamical complica- 
tions are hidden in the non-trivial Poisson brackets (5.7). 



VI. THE THERMAL PARTITION FUNCTION 

We are now in position to write down the classical par- 
tition function and compute (generally time-dependent) 



thermal expectation values. The classical p hase -space 
is defined by the initial conditions to Eqs. (5^4), and 
the ther mal weight is provided by the Hamiltonian in 
Eq. (|5.6|) ■ The correlation functions of the m agn etic fields 
A % a are then generated by (compare to Eq. (2.7)): 



V£fVA a l vw a s(g a ) 



x exp 



m + / d i xJf{x)A a l {x)\, (6.1) 



where A l a {x) is the solution to Eqs. ( J5.4D with the initial 
conditions {£f,A?,W a } (that is, Ef (t ,x) = £f (x), etc., 
with arbitrary to), and Q a a nd H , are e xpr essed in terms 
of the initial fields (cf. Eqs. (|J) and ft|)). 

It can be verified pa] t hat the phase-space measure 
T>£fT>AfVW a in Eq. |y) is invariant under the time 
evolution described by Eqs. (5.4), so that Z c i[J] is inde- 
pendent of the (arbitrary) initial time to, as it should. 
(This point is not trivial because of the non-canonical 
structure of the equations of motion.) 

Furthermore, strictly speaking, Eq. ( |6.l| ) is not well- 
defined without an UV cutoff (recall the discussion in 
Sees. II and III). This will be discussed in the next two 
sections. But be fore doing that, it might be illuminating 
to consider Eq. ( |6.l|) in some particular cases. 

For instance, for J = 0, it yields the same thermody- 
namic potential as the leading-order dimensional reduc- 
tion Specifically (compare to Eq. ( p,14| )): 



Z cl = JvA^VA* 



(6.2) 



where the fields Aq have been introduced as Lagrange 
multipliers, and the (Gaussian) functional integrals over 
£ f and W a have been explicitly performed. In particular, 
the integral over the auxiliary fields W a has correctly 
generated the Debye screening mass for the electric fields. 

Consider also the Abelian limit, where the classical 
theory can be exactly solved. Physically, this describes 
the propagation of soft electromagnetic waves in a hot 
QED plasma, with the fermions integrated over (in the 
HTL approximation). In the present approximation, the 
Abelian theory is quadratic, and does not require any 
UV cutoff. This is consistent with the fact that, in QED, 
soft and hard degrees of freedom correspond to different 
kind of fields (the hard fields are fermions, while the soft 
ones are photons), so, unlike QCD, there is no danger of 
overcounting. 

A straightforward, if lengthy, calculation yields (see 
the Appendix for details): 



Z c i [J%] 



exp 



d^yJiWDfUx-ypM} (6-3) 



with (compare to Eq. (E^q)): 
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DfAx-y) 



d 4 fc 
(2^ 



-ik-(x-y) 



■pij(k). 



(6.4) 



Here, Pij(k) is the photon spectral density in the HTL 
approximation and in the temporal gauge: 



Pij(k ,k) = (Sij - k i k J )p T (k ,k) + 



^p L (k ,k), (6.5) 
k 



with the (transverse and longitudinal) spectral functions 
Pt.l defined as standard They include both delta 
functions for on-shell quasiparticles and off-shell pieces 
corre spo nding to Landau damping. For instance (cf. 
Eq. Q): 



p T (fco , k) = 2 Im A T (k + ir), k) 

= 2Tre(k )z T (k)S(k 2 - Lu 2 (k)) + T (ko, k), (6.6) 

where zr(k) denotes the residue of the (time-like) poles at 
fco = ±^r(fc), and f3r(ko, k) oc ImHr(fco, k) has support 
only at space-like momenta, k^ < k 2 . 

More generally, in QCD, D cl will enter as a propagator 
in the perturbative expansion of Eq. (3.1). 



VII. THE HARD-SOFT MATCHING 



As already mentioned, Eq. (S.l) must be supplemented 
with an UV cutoff A which is meant to separate be- 
tween hard and soft gluons: gT «C A <C T. Corre- 
spondingly, the effective Hamiltonian must include some 
A-dependent counterterms, chosen so as to cancel the de- 
pendence on A in the calculation of physical quantities 
(a cancellation to be referred to as "matching"). The 
matching can be done only to one-loop order (since the 
HTL theory is already an one-loop approximation). But 
this is enough to remove the most dangerous, linear, UV 
divergences of the classical theory (see below) . 

For more clarity, consider first a very simple cutoff pro- 
cedure: the UV cutoff A is implemented as an upper 
momentum cutoff on the loop integrals of the classical 
theory^. (This would be just fine if we were interested 
in perturbative calculations.) For matching purposes, we 
need the A-dependent corrections to the soft amplitudes, 
to one-loop order in the effective theory. Such corrections 
are easy to isolate Jig] : they are precisely the "HTL's" of 
the classical theory. Namely, they differ from the usual 
HTL's onl y by the replacement of the correct Debye mass 
mjj, Eq. ( [2.12] ), wi th th e A-dependent classical "Debye 
mass" m%(A), Eq. fl2~ilj ). 

In this case, the matching is very easy: just replace the 
physical Debye mass m 2 D in the definition of the effective 



4 By "loops", I mean here, of course, the classical, thermal, 
loops, as generated in Eq. (3.1) by the non-linearities of the 
classical theory. 



theory (that is, in all the equations in Sees. V and VI) 
with the following A-dependent parameter: 



m 2 H (A) = m 2 D -m 2 cl {A) = 



g 2 C A T 2 



1 



Then, to one-loop order, the effective theory generates 
the correct soft amplitudes, that is, the usual HTL's with 
total strength m 2 H (A) + m 2 cl {A) — m 2 D . 

For non-perturbative calculations, however, one can- 
not use a simple momentum cutoff, as above, but one 
rather needs an UV cutoff which can be also implemented 
on the lattice. This cannot be the lattice spacing it- 
self: indeed, a finite (and relatively large) lattice spacing 
a ~ 1/A ^ 1/T will introduce lattice artifacts (essen- 
tially, the loss of rotational symmetry) which will make 
impossible the matching with the hard sector p3| , p4"| |. 

A first solution to this problem have been proposed in 
Ref. there, a smooth cutoff has been introduced in 
the classical theory by replacing, in Eq. (|5.6|), 



Tr WW - Tr W J [ ^ ) W 



(7.2) 



where f(z) = 1 + z 2 and the trace refers to color in- 
dices. In principle, such a higher-derivative cutoff can be 
indeed implemented on the lattice, as discussed in Ref. 
p4j . Then, the matching can be done (at least for the 
non-perturbative correlation functions) by adding a A- 
dependent counterterm to the Debye mass, as in Eq. (7.1) 
]18[ . In practice, however, this might be quite tedious, 
as it requires highly improved lattice Hamiltonians p4| . 
Fortunately, there exists a simpler way to perform clas- 
sical non-perturbative calculations without running into 
ultraviolet problems. This will be presented in the next 
section. 



VIII. THE EASY WAY 

Note first that, within the effective theory, it is easy to 
distinguish between the gluon fields which are responsible 
for the non-perturbative IR phenomena and those which 
are responsible for the leading UV behaviour: 

The non-perturbative phenomena are generated by soft 
and quasistatic magnetic gluons (cf. Sec. Ill), which 
are described by the off -shell piec e oc [3t of the gluon 
propagator (cf. Eqs. (6.4) and (3.£)). 

The UV behaviour, on the other hand, is determined 
by the relatively hard (k ~ A) magnetic gluons, which are 
essentially on-shell (and, actually, even on the tree-level 
mass shell) . Indeed, for momenta k ~ A ;§> run , the HTL 
corrections are negligible and the classical propagator in 
Eq. ( |6.4p reduces to its tree-level expression, Eq. ([Q|). 
In fact, the "classical HTL's" mentioned in Sec. VII (in 
relation with the matching problem) have been computed 
[ p3| , p~8 J24 1 by using the free propagator Df}(k) for the 
hard (k ~ A) gluons around the loop. 
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This suggests that one can avoid the UV problems of 
the classical theory by subtracting (in a way to b e spec- 
ified shortly) the tree-level piece po(k), Eq. fl2.8|), from 



the magnetic spectral density pT{k), Eq. (6.6). 

To understand how this can be achieved in practice, 
recall from Sec. VI that the magnetic two-point func- 
tion D^(k) = (T/ko)pT(k) is obtained by averaging the 
classical solution over the initial conditions {A l ,£ l ,W} 
at the (arbitrary) initial time to- By using the equations 
in the Appendix, it can be verified that, as to —oo, 
(i) the off-shell piece of the propagator, proportional to 
/?Tj is generated by the integral over the colour fluctu- 
ations W Q , while (ii) the on-shell piece arises from the 
integration over the initial colour fields A l a and £ % a . 

This is physically intuitive: (i) the off-shell piece de- 
scribes Landau damping, so it should arise from the cou- 
pling of the soft fields to the long-range colour fluctu- 
ations W a of the hard particles; (ii) at high momenta 
k ~ A ^> 771 jj , the one-shell piece goes into the tree- level 
propagator D^ik), which i s th e sam e a s in the purely 
Yang-Mills theory (cf. Eqs. (|J) and (U)), and is tnere ~ 
fore independent of the W's. 

This can be also verified by a rapid calculation (see the 
Appendix for more details): To this aim, it is enough to 
consider the linearized (or Abelian) version of the effec- 
tive theory, since this is what determines the propaga- 
tor. The corresponding equation for th e transverse field 



Ai r (x) reads, for x > t (cf. Eq. (|A.16| )) 



01 



V 2 )At+ r<i 4 yIlT(a : -y)Ai r (y)=&(x), (8.1) 

Jt 



where £ l is det erm ined by the initial fluctuations 
W(x,v), cf. Eq. (P 



= 771 



I I ^^ W (x-v(x -to),v) ) (8.2) 



oo, Eq. (8.1) can be easily solved by Fourier 



If t - 

transform. For vanishing initial conditions, A l T — £ ^ — 
0, the solution reads: 



Ai r (k) = A T (k)&(k), 



(8.3) 



wher e A x(k) is the retarded magnetic propagator in 
Eq. (3.1). The correlation function: 



(A^A^ip)) = (2Tr) 4 8^(p + k)D^(k) (8.4) 

can then be obtained by performing the functional aver- 
age o ver t he W' s. E quivalently, note that, according to 
Eqs. ( |6.lD and (5.6), the W's can be seen as Gaussian 
random variables with the local correlation function: 

(W(x,v)W(y,v')} = (T/m 2 D )5^( X -y)5(v,V). (8.5) 

This immediately implies: 

(e(x)?(v))=m 2 D T J ^W'Cx-y-vOro-tfo)), 

(8.6) 



and therefore, after some simple algebra, 

D#(u, k) = -2 (T/lu) \A T (k)\ 2 Imn T (w, k) 
= {T/u)(3 T {u,k), 



•7) 



where in writing the last equation I have recognized the 
imaginary part of the HTL polarization tensor |l[ : 



/dfl 
— v 1 v 3 5(lu - v • k 



together with the defin ition ( 3J: ) of the off-shell spectral 
function (3t- Eq. ( |3.7| ) is precisely the La nda u dam ping 
piece of the magnetic propagator in Eqs. (6.4)-(|6~6|). As 
anticipated, this has been generated here by averaging 
over the initial fluctuations W, while keeping constant 
the initial fields A l T and £ % T . 

The above discussion suggests the following strategy 
for non-perturbativc calculations in the classical effective 
theory: The (fully non-linear) equations of motion (5.4) 
should be solved with the following initial conditions at 
t = t Q : 

A a (x) = 0, £ a T {x) = 0, W a (x,v), 

ik£l{k) =m%J^W a (k,v). (8.9) 

That is, the initial colour fluctuations W a are arbitrary, 
but all the initial gauge fields are taken to be zero, except 
for the longitudinal co mpon ent of the electric field which 
is fixed by Gauss' law (5.5). The system is let to evolve 
for a long time (so as to simulate the limit to — » — oo), and 
then the correlation function of interest (e.g., Eq. (£J^)) 
is evaluated along the classical solution. This yields the 
correlation function as a functional of the W's, say 51 [W], 
which is then averaged as follows (below, Af = 1/(1)): 



(si)=N J v£ L vws{g)n[w] 



x exp 



£»*)}. (8,0) 



Clearly, this equation involves only a r estricted phase- 
space integration (compare to Eq. ( |6.l| )), in the sense 
that the initial transverse fields are kept constant (and 
equal to zero). But this is precisely what we need in 
order to generate the correct IR physics, without being 
sensitive to the hard classical modes. 

M ore precisely, the perturbative expansion of 
Eq. ( |S.10 ) will involve the reduced transverse propagator 
, Eq. ( |S.7D , rather than the full classical propagator in 
Eq. (|6.4|). The one- loop classical amplitudes constructed 



with will be UV finite, since dominated by soft loop 
momenta k ^ gT. For instance, the tadpole diagram 
constructed from Dj, reads: 



n 



d 4 fc 



D#(k) = T 



= T 



(2tt)4 

d 3 fc / 1 
(2tt) 3 [I? 



d 3 k 



2lTLU 



z T (k) 



.11) 
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which is indeed finite (and of 0(m£>T)) since the inte- 
grand in the last equation behaves like m 2 D /k A at mo- 
menta k gT j[o|. ( To ob tain this last equation, I 
have used the sum rule ( A. 21 ) for pr/u, together with 
Eq. (3.6).) This should be contrasted with the usual 
(classical) tadpole, as generated by the full propagator 
which is linearly UV divergent: 



n 



d 4 fc 

(2tt)* 



D$(k) = T 



d 3 fc 



T 



d 3 k 



1 



AT. 



.12) 



In general, all the "classical HTL's" of Sec. VII will be 
identically zero when computed with the reduced propa- 
gator Dy, so there is no need for matching. 

On the other hand, the reduced propagator generates 
the correct infrared structure, as determined by the off- 
shell, soft (fc ~ g 2 T) and quasistatic (ui ?5 g 4 T) mag- 
netic gluons, In fact, the spectral weight of the non- 
perturbative gluons is concentrated at low frequencies 
u <C k, so it is entirely contained in D^- 
call the low frequency behaviour of (3t / w (cf. 



To see th is, r e- 



Eq. 



P T (0J < k) 7T 



D ' 



2 /fc 6 + (TTm 2 D uj/4:) 2 



(U)): 
(8.13) 



For k <C rri£), this is strongly peaked at lu = 0, with a 
width Aw ~ fc 3 /m 2 D . (See Fig. [I] for an illustration.) For 
k ~ g 2 T, this yields u> ^ g 4 T. In particular, as k — » 0, 
@t{w, k)/ui — > (27r/fc 2 )J(tj), which is just another way to 
see the "dimensional reduction" |lO(l. 



14000 
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FIG. 1. The quantity pT(co,k)/ui as a function of ui for 
k = O.Iiiid. All the quantities are made adimensional by 
multiplying them by appropriate powers of mn- 



IX. CONCLUSIONS 

I have presented here a classical field theory for 
hot QCD which describes the non-perturbative infrared 



physics to lowest order in g and is well suited for numer- 
ical simulations on a lattice. This is an effective theory 
for the soft gluons (with momenta k gT), as obtained 
after integrating out the hard ones (fc ~ T) in the HTL 
approximation (i.e., to leading order in g). The classical 
theory has a built-in UV cutoff A (with gT < A < T) 
which can be introduced as a higher-derivative cutoff for 
lattice purposes. In practice, there is a simpler way to 
avoid UV problems, which is to restrict the phase-space 
averaging to the colour fluctuations of the hard gluons 
(the "auxiliary fields" W a (x,v)). 

The classical theory allows for three-dimensional lat- 
tice simulations of the real-time dynamics. In this sense, 
it appears as a generalization of the (leading-order) di- 
mensional reduction (PR ), to which it reduces in the 
static limit (cf. Eq. (3.2)). However — and unlike DR 
which, at least in principle, can be pushed to arbitrarily 
high orders [jl6{l7| — , the present procedure is not sys- 
tematic. That is, one cannot improve the classical the- 
ory by integrating out the hard modes, e.g., to next-to- 
leading order in g, not only because the latter operation 
would be technically involved, but also because such an 
improvement would not be consistent with the validity of 
the classical approximation. Different to say, the quan- 
tum corrections to the dynamics of the soft fields will 
generally enter on the same footing as the hard-field cor- 
rections to the HTL's. Besides, the whole construction 
of the classical theory as given above has relied in a very 
non-trivial way on the specific structure and symmetries 
of the HTL's. 

Still, the effective theory formulated as above is accu- 
rate enough to allow for leading-order calculations (e.g, 
on a lattice) of non-perturbative real-time correlations 
like the hot baryon number violation rate |13 LlJl9|,01, 
the quasiparticle damping rates [ 
coefficients (e.g., the bulk viscosity 

The effective theory should be further compared to the 
method developed in Refs. ]22| , where classical coloured 
particles have been introduced to simulate the HTL's on 
the lattice (see also Refs. [Q,^5[). The present formal- 
ism has the advantage to avoid the UV problem of the 
classical theory, thus allowing for the continuum limit 
a — > to be eventually taken. Moreover, this formal- 
ism should also include, as a special limit, the effective 
theory recently proposed by Bodeker ^?J. Therefore, by 
performing accurate lattice calculations in the present 
formalism, it should be possible to improve over the cal- 
culation in Ref . [^2| , and eventually distinguish between 
the various theoretical predictions ^3j,^4|j2^] for the hot 
baryon number violation rate. 



Appendix A 




or some transport 



Below, I shall derive Eqs. (6.3)-(6.4) for the thermal 
partition function in QED. The first step is to solve the 
Abelian version of Eqs. (5.4), namely (E l = —A 1 ): 
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{dl - V 2 )S ij + d l 8 3 A j (x) = m 2 D I W(x, v), 

47T 



(fib +v - V)W(z,v) = vE(ic), 

with the initial conditions at £q = 0: 

A^O.x) =A(x), i J (0,x) = -r(x), 
W(0,x,v) = W(x,v), 

which are constrained by Gauss' law: 



V ■£ 



, L) i ^W(x,v) = 0. 



(A.l) 



(A.2) 



(A.3) 



To completely fix the gauge, I further restrict the initial 
fields to be tranverse: V • A = 0. 

We thus have a linear initial value problem which can 
be solved by Laplace transform, in the standard way. 
(See, e.g., Ref. pj| for a similar problem.) 

Conside r fir st the Vlasov equation (i.e., the second 
equation (A.l)) for W(x,v). Its solution is conveniently 
written as: 



W = Wi, 



Wfi, 



(A.4) 



where Wind is the piece induced by the gauge fields, i.e., 
the solution to the Vlasov equation with zero initial con- 
dition (Wi n d(0, x, v) = 0), while Wfi satisfies the homo- 
geneous equation: 



(do + v ■ V) W fj = 0, 



(A.5) 



with the initial condition W/z(0, x, v) = W(x, v). It im- 
mediately follows that: 



W f i(x,v) = W(x- V£C ,V), 



(A.6) 



which is referred to as "the fluctuating piece" since it is 
determined by the initial charge fluctuations W. Also, 
for xo > t : 

W ind {x, v) = J d 4 y 0(yo) G R (x, y\v) v • E(y), (A.7) 

where Gr(x, y\v) is the retarded Green's function of the 
line derivative do + v • V : 

G R (x, y|v) - 9(x Q - vo) 5^ (x - y - vfo - y Q )) . (A.8) 



Eq. ( A.4 ) implies a similar decomposition for the current 
(5.3): f = j\nd + where is the fluctuating piece: 



/" dfi 

1 / ^u l W(x- va;o,v), 



(A.9) 



and j\ nd is the current determined by Wj 



ind ■ 



f dO . 

dV(yo)IF(a;-yL4%) 
*£> / ^uV\A J '(x-va5o), 



(A.10) 



where in writing the second line I have used Eq. ( |A. 7j ) to- 
gether with an integration by parts, and I have identified 
the (retarded) polarization tensor in the HTL approxi- 
mation: 



n" (./•-//') = nr v I ^± v i v id G R (x,y\v) 



d*k 



(2^) 4 

with the familiar HTL expression [jlj: 

dn vV 



IF (A) = m^fco / -r- 



47r /cq — v • k + i?7 



(A.ll) 



(A.12) 



The Maxwell equation (i.e., the first equation (A.l)) 
reads then (for xq > to): 



(d 2 - V 2 )^ + A?(x) = jt nd (x) + e(x), (A.13) 



and can be most easily solved by going to the three- 
momentum space and then separating the Fourier modes 
into transverse and longitudinal components: 



A*(t,x) = J 



d 3 fc 



A l (t,k) = A l T (t,\t) + k l A L (t,k), 



(A.14) 



where k ■ Ay = 0. A similar decomposition holds for the 
initial fields A 1 and £ l , with Al = (ga uge f ixing) and 
£l as determined by Gauss' law (cf. Eq. (A.3)): 



ik£ L (k) = m 2 D { _W(k,v). (A.15) 

Let me c onside r the transverse sector in more detail. 
Eqs. (1A.13D and ( |A.10| ) yield: 



(flg + fc 2 )A^(i) + / dt' 6{t')U T (t - t')A l T {t') 
-a\t) = &(t), 



(A.16) 



where the k-dependence is impl icit, £^(i, k) is the trans- 
verse projection of £*, Eq. (|A.9|), and 



a'(t, k) 



, /" ^e-^wV^Oc), (A.17) 



which is transverse as well: k • a = (N.B. A 1 = A l T .) 

The solut ion to Eq. (A.16) correspondin g to the initial 
conditions ( |A.2[) can be written as in Eq. ( kA ) : 



A r (t,k) = A? ld (i,k) + A^(t,k), (A.18) 

where A™ d is the piece determined by the initial electro- 
magnetic fields {A l T ,£^}, while Ajf is generated by the 
Ws (via £^). Specifically: 



A^ d (i,k) 



/dw , , r, 9 , coswi „ . i , . 
— p T {uj)\k 2 AT £ T sinujt} (A.19) 







and, respectively, 



A£(t,k) = J dt' 6(t') A T (t - 1') £ T (t') 



(A.20) 



wher e A y is the (retarded) transverse pro paga tor (cf. 
Eq. (3.1 )), wi th sp ectra l density px (cf. Eq. (|6.6|)). 

Eqs. ( A. 18 ) to ( A.20 ) express the transverse solution 
A-Tit, k) a s an explicit func tional of the initial condi- 
tions (A. 2) (recall Eq. ( A.9)). The fluctuating piece in 



Eq. ( A.2C ) is the solution to the inhomogeneous equa- 
tion ( A.16D with zero i nitial conditions: A l T = £ T = 0. 
The induced piece, Eq. (|A.19 ), satisfies the homogeneous 
equation (i.e., Eq. ( |A.lq ) with £ T = 0) with the proper 
initial conditions {A 1 ^-, To verify the latter, use the 
following sum-rules : 

f ^ ! \ 1 ;2 f dLJ g^M /A ni\ 

/ tt^Pt^) = 1 = k - . (A.21) 

J Zir J Zir oj 

The next step is to perform the average over the ini- 
tial conditions, a s ex pressed by the following functional 
integral (cf. Eq. (O)): 



Z cl [3 T ] = J V£ l VA l VW 5{G) e~ fm+ f d4xjT Ar 

= Z md [J T ] x Z fl [J T ], (A.22) 
where (recall that Al = 0): 

Z md [J T ] = J V£ T VA T e~ mT+ J d4lJrA?,i , 

Z fl [3 T } = Jv£ L VW5{g)e- fmL+ f dix3TA T 1 (A.23) 



and I have denoted 
1 



TLt 



d 3 x(s 7 



B' 



(A.24) 



with H — Ht + T~Cl- 

The first inte gral i n Eq. ( A.23| ) is straightforward and 
yields (cf. Eq. (|A~T9|)): 



z^[jT]=ex P |-y y__j T ( W ).j 5 



1 + — j)p T (u,k)p T (uj',k)}. (A.25) 



fi 

The second integral involves which is proportional 



to W. Specifically, Eqs. QA.9|) and (A.20) imply: 



— v T W(k,v) 



x Im 



(w,fc) 
CJ — v • k + irj 



By using this expression, together with the identity: 



exp 



T 
2m 2 D 



d 3 fc 

(27T) 5 



5(v,v')- 



47T 



47T 



£(k,v) 



/C*(k,v')k (A.27) 



one eventually obtains: 

z cl [j T ]Er d [j T ]x^[j T ] 



cxpi- 



d 4 fc p T (k ,k) 



(2tt)4 fc 



J T (fc)| 2 . (A.28) 



As expected, this provides the transverse piece of the 
classical two-point function in Eqs. (|6.^)- ( |G~5| ) . 
The corresponding longitudinal piece follows from: 

Z cl [J L ] = J V£ L VW5(g)e-> 3nL+ J dixjLAL , (A.29) 

where Al = A^ is the solution of the longitudinal equa- 
tion of motion, and is a purely fluctuating field: that is, 
it is fully determined by the initial charge fluct uatio ns W 
(since Al = and £l ix fixed by Gauss' law ( A.15| )). A 
calculation similar to the one above eventually yields: 



Z c i [Jl 



exp 



T 



d 4 fc p L {k) k 2 



(2^) 4 k k 2 



J L {k)\ 2 \. (A.30) 



Note that the would-be pole at ko — » (as typical for 
the temporal gauge) is actually innocuous in the classical 
theory, where the current conservation d^J^ = ensures 
that the longitudinal current also vanishes as ko — > 0: 
JL(k 0l k) = (k /k) J Q (k , k). 



Together, Eqs. (A.28) and (A.3C) provide the Abelian 
partition function presented in Eqs. (|6.3|)~(6.4). 
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